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Abstract. Recent use of Type la supernovae to measure acceleration of the universe 
has motivated questions regarding their optimal use to constrain cosmological param- 
eters Qm, 'a and wq. In this work we address the question: what is the optimal 
distribution of supernovae in redshift in order to best constrain the cosmological pa- 
rameters? The solution to this problem is not only of theoretical interest, but can be 
useful in planning supernova searches. Using the Fisher matrix formalism we show that 
the error ellipsoid corresponding to N parameters (for N — 1,2,3) has the smallest 
volume if the supernovae are located at N discrete redshifts, with equal number of su- 
pernovae at each redshift and with one redshift always being the maximum one probed. 
Including marginalization over the "nuisance parameter" Ai changes this result only 
trivially. 

1 Introduction 

Recent use of Type la supernovae (SNe la) as standard candles ]y|^||[lj pro- 
vided an opportunity to obtain the distance-redshift relation - and measure the 
cosmological parameters - without recourse to the "distance ladder". SNe la 
provide strong evidence for the acceleration of the universe, thus implying that 
a component called "dark energy" with strongly negative pressure (p < —1/3/?) 
dominates the energy-density of the universe. The most likely candidate for the 
dark energy is the vacuum energy, with p = —p. SNe la provide the luminosity 
distance - redshift relation and thus effectively determine the parameters £Im 
and (2a, energy densities in matter and vacuum energy scaled to the critical 
energy density. The best-fit value using the current dataset of ~ 50 SNe la and 
assuming a flat universe is Qm = 1 — — 0.3, with X-a error of order 0.1. 

Given the importance of SNe la as standard candles — namely, determin- 
ing the contents of the universe and probing the dark energy — it is important 
to consider what can be done to improve upon the constraints on cosmologi- 
cal parameters. Two obvious improvements would be increasing the size of the 
supernova sample and better control over the systematic errors. Going to red- 
shifts beyond one is also very high on supernova cosmologists' wish list because 
that way one can differentiate between the effects of dust and dark energy. It is 
believed that, barring unusual scenarios, dust would make faraway objects pro- 
gressively fainter, while the dark energy becomes inoperative at z > 1 because 
the universe is essentially matter-dominated at those redshifts. Finally, an im- 
portant requirement would be to have supernovae located throughout the probed 
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redshift range, in order to make sure that the magnitude-redshift relation traced 
out corresponds to the dark energy, and not evolution or dust. The planned 
supernova space telescope, SNAP would satisfy all of the above-mentioned 
requirements. 

This work considers supernova search strategies for the most accurate deter- 
mination of cosmological parameters Qm and Qa (and possibly the equation of 
state ratio of the dark energy, wq, where 'Q' stands for 'quintessence' @). To 
this end, we ask: given the cosmological parameters we want to determine, what 
is the optimal distribution of supernovae in redshift in order to best constrain 
those parameters? At first glance this problem may appear of purely academic 
interest since we are not free to put supernovae where we please. However, super- 
nova observers have considerable freedom in choosing redshift ranges for their 
searches, by using filters sensitive to wavelengths corresponding to spectra at 
observed redshifts. Moreover, the increased difficulty in observing high-redshift 
supernovae means that, even with great improvement in supernova detection 
and follow-up techniques, it can be as time-consuming to observe one supernova 
at, say, z — 1.4 as many z < 1 supernovae. Hence, an observer will have to 
decide how much telescope time is to be allocated to specific redshift ranges to 
best constrain the cosmological parameters. 

In this work we make three assumptions: 

• Magnitude uncertainty, a m , is the same for each supernova irrespective of 
redshift (this is a pretty good approximation, at least for current data sets). 

• Total number of observed supernovae is fixed (rather than the total ob- 
serving time, for example). 

• There is an unlimited number of supernovae at each redshift. 

None of these assumptions is required to use the formalism we present. More- 
over, the results we present should qualitatively not be very different from those 
obtained when the constraints above are relaxed. We make the assumptions 
above to illustrate our approach and simplify the analysis. 

2 The Most Accurate Parameter Determination 

We tackle the following problem: given n supernovae and their corresponding 
uncertainties, what distribution of these supernovae in redshift would enable 
the most accurate determination of cosmological parameters? In case of more 
than one parameter, we need to define what we mean by most accurate determi- 
nation of all parameters simultaneously. Since the uncertainty in measuring N 
parameters simultaneously is described by an iV-dimensional ellipsoid (at least 
under assumption that the total likelihood function is gaussian), we make a 
fairly obvious and, as it turns out, mathematically tractable requirement that 
the ellipsoid have minimal volume. This requirement corresponds to the best 
local determination of the parameters. 

We now show that volume of the ellipsoid is given by 



l/ocdet^)- 1 / 2 , 



(1) 
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where the sign of proportionality expresses our ignorance of a numerical factor, 
and F is the Fisher matrix [||f7| 

where L is the likelihood of observing data set y given the parameters p\ . . .pjy. 
Although this relation might be familiar/obvious to mathematically inclined 
cosmologists, we present its derivation for completeness. 

To prove equation (|l|), consider a general uncertainty ellipsoid in n-dimensional 
parameter space. The equation of this ellipsoid is 

X T FX = 1, (3) 

where X = {x\Xi . . . x n) is the vector of coordinates and F the Fisher matrix. 
Let us now rotate the ellipsoid so that all of its axes are parallel to coordinate 
axes. Equivalently we can rotate the coordinates to achieve the same effect, by 
writing X rot = UX, where U is the orthogonal matrix corresponding to this 
rotation. The equation of the ellipsoid in the new coordinate system is 

X£ t FX Iot = l, (4) 

or equivalently, in the original coordinate system 

X T F mt X = l, (5) 

where F rot = U T FU is the Fisher matrix for the rotated ellipsoid, and has the 
form F pal = diag(l/al, ... with <n the i-th axis of the ellipsoid. The 

volume of the (rotated) ellipsoid is obviously 

N 

V ocn^^det^ot)- 1 / 2 . (6) 

Then, since det(F) = det(F rot ) and rotations preserve volumes, we have 

V oc det(F mt )^ 2 = det(F)- 1 / 2 , (7) 
and this completes the proof. 



2.1 Fisher Matrix for Supernovae 

To minimize the volume of the ellipsoid we therefore need to maximize det(-F). 
Fisher matrix for the case of supernova measurements was first worked out by 
Tegmark et al. || , and we briefly recapitulate their results, with slightly different 
notation and one addition. The measurements are given as 

m„ = 5 log [H d L (z n , f2 M , Oa)] + M + e n , (8) 
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where m n is apparent magnitude of the nth supernova in the sample, cIl is its 
luminosity distance, M. = M — 51og(i?o) + 25 (with M absolute magnitude 
of a supernova), and e„ is the error in that measurement (assumed to be drawn 
from a gaussian distribution with zero mean and standard deviation <x m ). Note 
that M contains all dependence on Hq, since <1l oc 1/Hq. 
The Fisher matrix is given by [pi 



1 N 



(9) 



where w's are weight functions given by 



Wi{z) 



kS'[kI(z)} 



dl I{z) 
dpi 



In 10 [ S[kI(z)] 
if the parameter pi is flu or fl A (or £2q), or else 
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kS'[kI(z)] dl 
S[kI{z)} Wi 



if pi is wq . Also 



iJ d L = (1 + 2) 



S(reJ) 



1 

2^2 



(10) 



(11) 



(12) 



fsmh(x), if ^tot > 1; 
x, if ^tot = 1; (13) 
sin(a;), if /?tot < 1- 

I(z,n M ,n A ) = f A(z')- 1/2 dz' (14) 
Jo 

A(z) =n M (i + z) 3 + n A + K 2 {i + z) 2 (15) 

K 2 = 1 - O m - (16) 

Here Qm-i and k 2 are the energy densities in matter, cosmological constant 
and curvature respectively divided by the critical density, and ^tot = f^M + ^A- 
Later on we also consider a more general equation of state for the exotic energy, 
and replace the second term in equation ( |l5| ) by fiq{l + z) 3 ( 1+WQ h 

In addition to flu-, and wq, the magnitude-redshift relation also includes 
the "nuisance parameter" Ai, which is a combination of the Hubble parameter 
and absolute magnitude of supernovae, and which has to be marginalized over in 
order to obtain constraints on parameters of interest. Ignoring A4 in the Fisher 
matrix formalism (that is, assuming that M. is known) leads to a serious un- 
derestimate of the uncertainties in other parameters (of course, fairly accurate 
knowledge of Hq could be used to obtain M from a sample of nearby supernovae, 
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and thus determine M.). We continue to ignore M. for simplicity and mathemat- 



ical clarity, and in section 3.4 we show that including marginalization over A4 
changes our results rather trivially. 

The Fisher matrix can further be written as 

AT /• 2max 

F%j = — / g(z)Wi(z)wj(z)dz, (17) 



where 



1 N 

= ( 18 ) 



is the (normalized) distribution of rcdshifts of the data and z max is the highest 
redshift probed in the survey. Our goal is to find g{z) such that det(F) is maxi- 
mal. g(z) is essentially a histogram of supernovae normalized to have unit area. 
Note that the maximization of det(F) will not depend on N and a m , so we drop 
them from further analysis. To consider a non-constant error o~ m (z), one simply 
absorbs a m (z) into the definition of weight functions w(z). 



3 Results 

3.1 Determination of One Parameter 

We first consider the trivial - but instructive - case of one parameter. Then we 
need to maximize J Zmax g(z) w\(z) dz, subject to J Q Zmax g(z) dz = 1 and g(z) > 0. 
It is quite obvious that the solution is a single delta function at the redshift where 
ui\(z) has a maximum. For any given parameter w\(z) will have a maximum at 
Zmax an d that is where we want all our supernovae to be. This result is hardly 
surprising: we have a one-parameter family of curves m(z), and the best way to 
distinguish between them is to have all measurements at the redshift where the 
curves differ the most, at z max . 

As an example, Fig. [j] shows magnitude-redshift curves for the fiducial Qm = 
0.3 model with the assumption Qa = 1 — f^M (flat universe). As Qm is varied, 
the biggest difference in m(z) is at the highest redshift probed. In order to best 
constrain Qm, therefore, all supernovae should ideally be located at z = 1.0, our 
assumed maximum redshift. 

3.2 Determination of Two Parameters 

A more interesting - and relevant - problem is minimizing the area of the el- 
lipse which describes the uncertainties for two parameters. The expression to 
maximize is then 

{z)w\(z)dz I g{z)u>2{z)dz — ( / g(z) W\(z) u>2{z)dz 

z=0 J z=0 \Jz=0 
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g{zi)g{z 2 )w (z 1 ,z 2 )dz 1 dz 2 , (19) 

'zi=0 Jz 2 =0 
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Fig. 1. Dependence of the magnitude- redshift curve on a single cosmological parameter. 
We assume a fiat universe and vary Qm = 1 — Qa- Magnitudes for the Qm =0.3 fiat 
model were subtracted from both curves for easier visualization 

where w(z\, Z2) = wi(zi)w2(z2) — w\ (22)^2(^1) is a known function of redshifts 
and cosmological parameters (see Fig. @) and g{z) is subject to the same con- 
straints as in case of one parameter. 
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Fig. 2. Left panel: Function w 2 (zi, 22) for the case when Qm = 0.3 and Qa = 0.7. Right 
panel: Variation of the magnitude-redshift curve when two cosmological parameters, 
Qm and Qa, are varied independently 



Despite a relatively harmless appearance of expression ( |l9| ) , we found it im- 
possible to maximize it analytically. Fortunately, for all practical purposes nu- 
merical solution will be sufficient. Returning to the discretized version of equa- 
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tion (fi7|), we divide the interval (0,z max ) into B bins with giN supernovae in 
bin i. Then we need to maximize 

B 

^2 g i g j w 2 {z i ,z j ) (20) 
i.i=i 

subject to 

B 

^<?* = 1 and 9l >0. (21) 

i=l 

Equations ( |20| ) and ( pl| ) define a quadratic programming problem — extrem- 
ization of a quadratic function subject to linear constraints. Since w 2 (zi,Z2) is 
neither concave nor convex (see Fig. ^|), we have to resort to brute force maxi- 
mization, and consider all possible values of gi. The result of this maximization 
is that the optimal distribution is two delta functions of equal magnitude: 

g(z) = 0.50 6(z - 0.43) + 0.50 6(z - 1.00), (22) 

where all constants are accurate to 0.01. Thus, half of the supernovae should 
be at the highest available redshift, while the other half at about 2/5 of the 
maximum redshift. 

This result is not very sensitive to the maximum redshift probed, or fiducial 
parameter values. If we increase the maximum available redshift to z max = 1.5, 
we find two delta functions of equal magnitude at z = 0.57 and z = 1.50. If 
we change the fiducial values of parameters to flu = 0.3 and Qa — (open 
universe), we find delta functions of equal magnitude at z = 0.47 and z = 1.00. 

Finally, let us consider a different choice for the two parameters — for ex- 
ample, f?M and wq (equation of state of the dark energy) , with fiducial values 
Qm = 0.3 and wq = — 1 and with the assumption of flat universe (J?q = 1 — J?m)- 
Then, assuming z max = 1.0, we get 

g(z) = 0.50 6(z - 0.36) + 0.50 S(z - 1.00), (23) 

and again optimal distribution of supernovae is similar to the case of flu and 
Qa as parameters. 



3.3 Determination of Three Parameters 

We further consider parameter determination with three parameters J?m, Qq 
and w. Elements of the 3x3 Fisher matrix are calculated according to expression 
(p7|), and we again maximize det(F) as described above. The result is 

g{z) = 0.33 5{z - 0.21) + 0.34 5{z - 0.64) + 

0.33 8{z - 1.00), (24) 

with all constants accurate to 0.01. Hence we have three delta functions of equal 
magnitude, with one of them at the highest available redshift. 
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It appears impossible to prove that TV parameters are best measured if the 
data form TV delta functions in redshift. However, it is easy to prove that, if 
the data do form TV delta functions, then those delta functions should be of 
equal magnitude and their locations should be at coordinates where the "total" 
weight function (e.g. w 2 (zi, z-i) in case of two parameters) has a global maximum. 
In practice, the relevant number of cosmological parameters to be determined 
from the SNe la data is between one and three, so considering more than three 
parameters is less relevant for practical purposes. 



3.4 Including Marginalization over A4 

So far we have been ignoring the parameter A4, assuming that it is known 
(equivalently, that the value of Hq and the absolute magnitude of supernovae 
are precisely known). This, of course, is not the case, and M. is marginalized 
over to obtain probabilities on cosmological parameters. Fortunately, when M. 
is properly included, our results change in a predictable and rather trivial way, 
as we now show. 
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Fig. 3. Magnitude-redshift curve when a single cosmological parameter, Qm, is varied. 
Marginalization over M is now included, which means that the curves are allowed to 
slide vertically as well 



Including Ai as an undetermined parameter, we now have an TV+l-dimensional 
ellipsoid (TV cosmological parameters plus M), and we want to minimize the vol- 
ume of its projection onto the TV-dimensional space of cosmological parameters. 
The equation of the TV-dimensional projection is 



X Fp T0 jX — 1 



(25) 
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and -Fproj is obtained as follows: 1) Invert the original F to obtain the covariance 
matrix F^ 1 2) pick the desired NxN subset of F _1 and call it F~ r ^, and 3) 
invert it to get F pl0 j. 

Minimizing the volume of the projected ellipsoid we obtain the result that 
the optimal supernova distribution is obtained with N delta functions in redshift 
obtained when ignoring M., plus the delta function at z = 0. All N + 1 delta 
functions have the same magnitude. The intuitive explanation for this result 
is illustrated in Fig. |], which shows the magnitude-redshift curves when Qm 
and M are varied (flat universe is assumed). This figure is the same as Fig. |l], 
except the curves are now allowed to slide vertically as well, corresponding to 
the variation in A4 . There are two locations of largest departure when the two 
parameters are varied, namely z = and z = 1. It makes sense then that those 
are precisely the locations where the supernovae should be, and our analysis says 
that we ideally need equal number of supernovae at each location. 

3.5 Optimal vs Uniform Distribution 

Are the advantages of the optimal distribution significant enough that one should 
consider them seriously? In our opinion the answer is affirmative, as we illus- 
trate in the left panel of Fig. |[ This figure shows that the area of the Qm-^A 
uncertainty ellipsoid is more than two times smaller if the SNe have the optimal 
distribution in redshift as opposed to the case of uniform distribution. 




Fig. 4. Left panel: Uniform (dark-blue) vs optimal (light-red) distribution in redshift. 
Shown are constraints on Qm and Qa when M was marginalized over. To get the 
absolute sizes of the ellipses, N = 100 SNe were assumed with individual uncertainties 
of a m = 0.15 mag. Right panel: Thinnest possible ellipse for given N and er m (dark- 
blue) is infinitely long in one direction. However, the smallest-area ellipse (light-red) is 
almost as thin, and for all purposes can serve as the thinnest ellipse 
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3.6 Requiring the Thinnest Ellipse 

If wc are using SNe la alone to determine the cosmological parameters, then we 
clearly want to minimize the area of the error ellipse (we consider the case of 
£2m and ft a as parameters in this section). However, supernova measurements 
will also be combined with other methods to determine the parameters. A prime 
example of a symbiosis between two or more methods is combining CMB mea- 
surements with those of supernovae (l^,|8| . These methods together can improve 
the determination of Qm and Qa up to a factor of 10 as compared to either 
method alone due to breaking of the degeneracy between the two parameters. 
As can be seen in Fig. 2 of ref. H, in combining the CMB with SNe la data 
one might hope for the thinnest ellipse possible coming from supernova measure- 
ments. Here by "thin" we mean that the combination Qm — is accurately 
determined. 

Finding the thinnest ellipse is the problem that we can solve using our for- 
malism, since the length of each axis of the ellipse is proportional to the inverse 
square root of an eigenvalue of the corresponding Fisher matrix. All we need to 
do then is maximize the larger eigenvalue of F with respect to the distribution 
of the supernovae g{z). 

The result is perhaps not at all surprising: to get the thinnest ellipse, all 
supernova measurements should be at the same (maximum) redshift, which at 
the same time implies an infinitely long ellipse. More generally, we find that 
changing the redshift distribution of supernovae doesn't change the thickness of 
the error ellipse greatly, but does change its length. By attempting to obtain a 
thinner ellipse, we end up only making it longer. Fortunately, the smallest area 
ellipse we found can in practice serve as the thinnest ellipse, as shown in the 
right panel of Fig. ||. 

3.7 Reconstruction of the Potential of Quintessence 

It has been shown recently that a sufficiently good supernova sample could be 
used to reconstruct the potential of quintessence out to z ~ 1 [ fl4]]l5| l . More gen- 
erally, equation of state ratio of the missing energy, wq, can also be reconstructed 

s 

In the spirit of our analyses above, we ask: what distribution of supernovae 
in redshift gives the smallest 95% confidence region for the reconstructed po- 
tential V((f>)7 To answer this question, we perform a Monte-Carlo simulation by 
using different distributions of supernovae and computing the average area of 
the confidence region corresponding to each of them. 

Uniform distribution of supernovae gives the best result among the several 
distributions we put to test. This is not surprising, because reconstruction of 
the potential consists in taking first and second derivatives of the distance- 
redshift curve, and the most accurate derivatives are obtained if the points are 
distributed uniformly. For comparison, gaussian distribution of supernovae with 
Z = 0.7, a z = 0.4 gives the area that is 10 — 20% larger. 
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4 Discussion 

We considered supernova search strategies that produce the tightest constraint 
on cosmological parameters by minimizing the volume of the error ellipsoid cor- 
responding to those parameters. We first proved that, assuming that the total 
likelihood function is gaussian, the volume of an iV-dimensional error ellipsoid is 
proportional to the inverse square root of the determinant of the corresponding 
Fisher matrix. Then we showed that, if the supernova measurements are used 
to determine N parameters with N = 1,2 or 3, this volume is minimized if 
the distribution of supernovae in redshift is given by N delta functions of equal 
magnitude. In particular, J?m and fl\ have the smallest error ellipse if half of 
the measurements are at z max , while the other half are at roughly 2/5z max . If 
M. is marginalized over, we need an additional supernova sample at z = (a 
total of N + 1 delta functions in redshift) . 

Our approach is quite flexible and can be applied in practice. For example, 
given the redshift-dependence of the measurement uncertainties for a given ex- 
periment, cr m (z), as well as the time it takes to obtain a measurement at a given 
redshift (determined by the telescope specifications), one can perform analysis 
similar to that in section ^ to infer the optimal search strategy. 

It is important to keep in mind that the best determination of the parameters 
around their fiducial values is only one of the possible objectives. As mentioned 
in the introduction, tracing out the magnitudc-rcdshift curve throughout the 
redshift range is important to test for systematics such as evolution or dust. A 
sample of nearby supernovae would be useful to further check for systematics 
(at least one low-z SN search program is already under way). Ultimately, the 
chosen supernova search strategy should combine all of these considerations. 
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